A New Approach to Quintessence and a Solution of Multiple Attractors 



Shuang-Yong Zhou^ 

^Department of Physics, Shandong University, Jinan, 250100, P. R. China 
Email : zhousyOmail . sdu . edu . cn 

We take a new approach to construct Quintessential models. With this approach, we first easily 
obtain a tracker solution that is different from those discovered before and straightforwardly find 
a solution of multiple attractors, i.e., a solution with more than one attractor for a given set of 
parameters. Then we propose a scenario of Quintessence where the field jumps out of the scaling 
attractor to the de-Sitter-like attractor, by introducing a field whose value changes a certain amount 
in a short time, leading to the current acceleration. We also calculate the change the field needs for 
a successful jump and suggest a possible mechanism that involves spontaneous symmetry breaking 
to realize the sudden change of the field value. 



Recent observations and experiments strongly indicate 
that the universe is spatially flat and currently undergo- 
ing accelerated expansion P, 0, Q . A negative pressure 
energy component, termed dark energy, is suggested to 
be responsible for the acceleration. The simplest candi- 
date for dark energy seems to be a positive cosmologi- 
cal constant, which is conventionally associated with the 
quantum vacuum energy. However, it is very tiny, com- 
pared with typical particle physics scales, which is the 
so-called fine-tuning problem [J]. It also suffers the so- 
called coincidence problem Q. Rather than dealing di- 
rectly with the dark energy a cosmological constant, var- 
ious alternative routes have been proposed, which usu- 
ally invoke dynamical scalar fields, such as Quintessence 
BMBBIS. Phantom ^ and Quintom [l3]. 

Quintessence invokes an evolving canonical scalar field 
slowly rolling down its potential (to some extent like the 
inflaton which drives the inflation in the early universe) 
with equation of state > — f . Motivated from obser- 
vational data [l^, [T3|, the Phantom invokes a negative 
kinetic energy with effective equation of state wph < — 1, 
having led to many interesting phenomena [T5| . 

Among the various Quintessential models, tracker so- 
lutions have attracted a lot of attention. The tracker 
field has an equation of motion with attractor-like so- 
lutions in the sense that a very wide range of initial 
conditions rapidly converge to a common, cosmic evolu- 
tionary track of Ptf,{t) and w^{t). The tracking behavior 
with < Wm occurs when T > 1 and is nearly con- 
stant (d(|r - l|)/dlna < |r - 1|), where F is defined 
as VV" /V'"^, with V the potential and ' the derivative 
w.r.t. the field [1^1 ■ It has been found that the general 
inverse power-law {V{(j)) — and exponential 

{V{(i)) = Vq exp(l/0)) potentials are tracker solutions (we 
have chosen = 8ttG =1). 

Another important class of Quintessential models are 
scaling solutions [13, [H, [l^, HO, IMI in which the energy 
density of the scalar field mimics the background fluid 
energy density. Namely scaling solutions are character- 
ized by the relation cx whose simplest realization 
is the exponential potential V^e'^'^ . As long as the scal- 
ing solution is the dynamical attractor, for any generic 
initial conditions, the field would sooner or later enter the 



scaling regime, being sub-dominant during radiation and 
matter dominated eras to satisfy the tight constraints 
from nucleosynthesis and structure formation, thereby 
opening up a new line of attack on the fine-tuning prob- 
lem 22]. However, exit from the scaling regime is needed 
so as to give rise to recent acceleration. 

The double exponential potential [1, [2^ of the form 

V{ct^)^V^{e-^^ + e-''*) , (1) 

provides a simple realization of the exit from the scaling 
regime. Such potentials can arise as a result of compact- 
ifications in superstring models. By properly choosing /i, 
V and initial conditions, one term in the potential dom- 
inated over the other before nucleosynthesis, giving rise 
to the scaling solution, while the situation has reversed 
recently, giving rise to a de-Sitter-like acceleration. How- 
ever, whether it is possible to obtain required values of /i 
and V remains a problem. In [9|], the authors considered 
the potential 

V[c^) = ^0 [coMpcj)) - 1]" , (2) 

which has two interesting asymptotical regions. One of 
these with {\p,(f)\ ^1, < 0) gives the scaling solution, 
while the other with (|/i(/)| <C 1), according to virial theo- 
rem, gives current acceleration with average equation of 
state {wif,) = [n — l)/(n -I- 1). As current data favor an 
equation of state close to —1 0, n should be close to 
0, which is mathematically viable, but seems unnatural 
physically. See [13, UHl for another two popular models. 

On the other hand, there is an attempt to search for a 
solution of two scaling regimes by coupling Quintessence 
to the matter [2^. Nonetheless, this scenario has faced 
severe challenges since it has been showed that it cannot 
be realized for a vast class of scalar field Lagrangians p3| ■ 

In this letter, we take a new approach to construct 
Quintessential models. With this approach, instead of 
proposing an interesting Quintessential potential directly, 
we first propose a relation between two quantities, F and 
A (defined as —V'/V), and then figure out the potential. 
First, we show that a tracker potential which is differ- 
ent from that discovered before can be easily obtained. 
Then we find it straightforward to get a solution of mul- 
tiple attractors, that is, a solution with more than one 
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attractor for a given set of parameters. In the particular 
case given in this letter, we have a scaling attractor and 
a de-Sitter-hke attractor. We thus propose a model in 
which the universe first evolves to the scaling attractor, 
and then, by introducing a field whose value changes a 
certain amount in a short time, the universe jumps out 
to the de-Sitter-like attractor to give the current acceler- 
ation. We also calculate the change the field needs for a 
successful jump and justify the introduction of this kind 
of field. 

To start, we consider the action of Quintessence (e ~ 1) 
(or Phantom (e = —1)) minimally coupled to gravity. 



S = 



-9\ 



(3) 



where we use the metric signature (—,+,+,+) and 
(V(/))^ = g^'^ d^4'du4>- In the flat Friedmann-Robertson- 
Walker spacetime, the equation of state for the 
Quintessential field (j) is given by 



P4, _ e</)2 - 21/(0) 



(4) 



P4> ^ 2V{(j)) 

The variation of the action ^ with respect to gives 

£0 + 3eH^ + y = . (5) 

Since we carry out cosmological dynamics of the 
Quintessential field (j) in the presence of a barotropic fluid 
whose equation of state is given by Wm — Pm/ Pm (in this 
paper, we assume that Wm is constant), Einstein equa- 
tions reduce to 



3 ^2' 



H=-\ [£0^ + (l+Wm)p„]. 

Introducing the following dimensionless variables 



r _ vv" 



(6) 
(7) 



V 



(8) 



Eq . ^ (111), ([7]) can be recast in the following form [13, 

da: o , ^ A 2 

diV = -3-+— ^^2/ 

+ ^a;[(l - w^)ex' + (1 + Wm){l - 2/')] , (9) 



dy _ _\/6 
dN ~ 2 

dA 
diV 



z 

+ ^y[(l - w,r,)ex' + (1 + m^)(l - y^)] , (10) 



(11) 



where iV = Ina (a is the scale factor), together with a 
constraint equation 



2 I 2 I Pm -. 

ex +y +3^ = 1 



(12) 



The equation of state and the fraction of the energy 
density flj, for the field (j) a-re, respectively, 



y 



(13) 
(14) 



To warm up, we note that for many Quintessential (or 
Phantom) potentials F can be written as a function of A. 
Let us take the case of Phantom potential of the form 



y(0) 



for example. It is found 



F = 1 



Vo 



[cosh{<T(j))Y 



1 

n 



na 



(15) 



(16) 



Substituting Eq. p6|) into Eq. pT|) , we can perform three- 
dimension dynamical analysis of the autonomous system. 
For a barotropic fluid background, there is a unique sta- 
ble fixed point (x = 0, y = 1, A = 0), which is a de- 
Sitter-like dominant attractor. For the case n = 1, it 
confirms the results of 0]. Note that we neglect the 
cases with ?/ < 0, as the system is symmetric under the 
reflection (x, y) — > (x, — y) and time reversal t — > —t. 

The direct way to get a Quintessential (or Phantom) 
model is to conceive (usually fairly carefully) a poten- 
tial that meets constraints from observations and exper- 
iments. However, encouraged by what has been showed 
above, let us take another route. 

We note that the dynamical system ^ [TOl fTT|) , is obvi- 
ously autonomous except for F. In fact, since the poten- 
tial V{4)) is only a function of the field </>, by the definition 
([5]), A and F can be written as 

A = P(0) , F = Q(0) . (17) 

If the inverse function of P{(t)) exists, then we have 

T^Q{p-\\))^f(\). (18) 

It is noteworthy that in principle we can figure out the 
potential as a function of the field 0. Using the definition 
of A and F, Eq. (fT8| can be can rewritten as 



(19) 



V' V 
V" = —J{--)^F{VX). 



Let = V, then we get 

dft, _ 1 
W ~h 



F{V,h). 



(20) 
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Having figured out h{V), we can solve V'{(f)) = h{V{(j))) 
to obtain the potential V{(j)). Thus we can perform three- 
dimension dynamical analysis of the system ([U [TOl fTTj) 
with a fairly large amount of potentials beyond the ex- 
ponential case where the dynamical system reduces to 
two-dimension autonomous system. 

In the viewpoint of F as a function of A and consider- 
ing the powerful theorem presented by (l6j . it is easy to 
obtain a tracker field. As an example, we write T as 

r = i + ^, (21) 

which can be solved to give the potential 

= Voe^^'^'^+^^ , (22) 

where Vq{> 0) and (3 are integral constants. We note 
that it is different from the general inverse power-law 
(1/(0) = Ca/(l>") or exponential (F(0) = Vq exp(l/0)) 
potentials. 

Obviously, F > 1 if a > 0. To confirm this is a 
real tracker solution, we perform the condition d(|F — 
l|)/dA^< |r- 1|, and we get 

(23) 

Substituting Eq. JTl]) for Eq. ([23|, we obtain 

|aa;|<|A|. (24) 

Considering the tracking condition |A| ^ \^/x\ we 
finally get 

(25) 

Tracking behaviors exist for a wide range parameters 
and initial conditions, which solves the fine-tuning prob- 
lem. However, due to the w-il relation that alleviates the 
coincidence problem, it is difficult to obtain current equa- 
tion of state Wo < —0.8. A numerical solution of the cos- 
mic dynamical evolution with tracking behavior, where 
for simplicity we have neglected the matter-dominated 
era, is given by Fig. [TJ 

Having witnessed the utility of this approach, we would 
like to go further. As showed below, we find it straight- 
forward to parameter F as function of A to get a so- 
lution of multiple attractors, i.e., a solution with more 
than one attractor for a given set of parameters. In the 
particular case given below, we have a scaling attrac- 
tor and a de-Sitter-like attractor and it is worth noting 
that two scaling solutions are problematic [27|. Thus 
we are encouraged to consider a scenario that the ini- 
tial conditions of the cosmic scalar field arc in the basin 
of a scaling solution and first the field evolves toward 
the scaling attractor. Then recently, the field jumps out 
to the basin of a de-Sitter-like dominant attractor, giv- 
ing rise to the current acceleration. In this scenario, the 
mechanism of exit from the scaling regime is different 
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FIG. 1: Evolution of W4, (red daslied line) and 57^ (green solid 
line) of r = 1 + a/X'^ {V{<j)) = Y^e°"^i<t'+I^)l'^ ^ a is chosen as 
2.8) with respect to A?^ = In a in the background fluid with 
Wm ~ 1/3. We choose initial conditions as Xi = 0.1, yi — 0.36 
and Xi — 17.8. For simplicity we have neglected the matter- 
dominated era. 

from those mentioned in the introduction, which typi- 
cally invoke fairly carefully conceived potentials with two 
asymptotical behaviors corresponding to the scaling case 
and the de-Sitter-like case respectively. Therefore attrac- 
tors in those models are not exact. On the contrary, the 
two attractors considered below are exact and we sug- 
gest some other physical reason to realize the exit from 
the scaling regime, rather than connect the two interest- 
ing regimes with more or less contrived potentials. The 
physical reason is formulated as the sudden change of the 
field value. 

Considering Eq. (|lip . we parameter F as 

r =1 + 1/(3+^ . (26) 
A 

There are at least the following two fixed points: 

• Point (a): 

{x = — \/6q/3, y = \/l — A — —a/3) is a 

de-Sitter-like dominant attractor, in which = 
— 1-1- /3. The eigenvalues of the Jacobi matrix 
of the dynamical system are 

Ail = 

Ai2 = -3H — , 

fi3 - -3(1 + Wrn) + a^/?^ • 

It exists if a^/?^ < 6 and is stable if < 3(1 + 
Wm) and P > 0. 

• Point (6): 

{x ^~V3{1 + Wrn)/V2ap, y =^3(l-0/2a2/32^ 
A = —a/3) is a scaling attractor, in which w^ = Wm 
and fi^ = 3(1 -I- Wm) / c^^ P"^ ■ The eigenvalues of the 
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FIG. 2: Evolution of x (red solid line), y (green dot-dashed 
line), A (blue dashed line) ofr = 1 + 1/ (3 + a/ X with respect 
to = In a in the background fluid with Wm = 0. a is 
chosen as —2.6, (3 chosen as 2. We choose initial conditions 
as Xi — 0.2, Ui = 0.36 and Xi = —1.3 for the thin lines, and 
Xi = 0, = 0.06 and Xi — —1.4 for the thick lines. Note that 
the attractor is a stable spiral. 



Jacobi matrix are 

3(1 + Wm) 



M2,3 = -|(1 - Wjn){l ± 



It exists ifa^/32 > 3(1- 
3(1 + Wjn) and /? > 0. 



-7-9w„, + 24(l + m^)2 
'm) and is stable if > 



For the stabilities of the fixed points a or b, we choose 
/3 > 0. We note that the fixed points a and b are typical 
for general scaling solutions and cannot be both stable 
for a given set of parameters [H, [M] . 

Besides, we find that the de-Sitter-like dominant fixed 
point (x = 0, y — I, A = 0) is stable, i.e., a de-Sitter- 
like dominant attractor which we denote as Point (c). It 
cannot be simply seen from the eigenvalues of the Jacobi 
matrix (/ii — — 3(H-w,„), ^2 ~ —3, /is = 0), since /ia = 

0. However, it can be seen from numerical simulation of 
the dynamical system, see Fig. [2] for a working example. 
It is found that when a < (a > 0), the region of A < 
(A > 0) in the phase space of the dynamical system is the 
basin of the attractor (x = 0, y = 1, A = 0), while the 
rest region is the basin of the attractor with A = —a/3, 

1. e., (a) or (6). This is desirable, since the basins are 
divided by a plane of equal A in the phase space, having 
nothing to do with x and y. Thus the initial values of x 
and y can be arbitrary. 

The corresponding potential of this case is 



Vo 



(27) 



where Vb(> 0) and are integral constants. For the 
stability of the potential, that is, the potential should be 



bounded below, we should choose [3 as 2,4,6,... and to 
obtain interesting cases we choose 77 < 0. 

Now we come to consider the scenario that the field 
exits the scaling regime to the de-Sitter-like regime due 
to a sudden change of the field value. To this end, we 
conceive 6 as 



with 



lit) 



^ = fit)^, 



1 t<t^ 
fj t>t 



1 ' 



(28) 



(29) 



where t is the cosmic time and fj is a constant with 
the subscript j representing some recent time when the 
field jumps. Note that f{t) is not necessarily of the form 
above, but it should have a certain amount of change of 
its value in a short time so that the field will not evolve 
back to the scaling attractor. We will first calculate the 
change the field needs for a successful jump and then sug- 
gest a possible mechanism to realize the sudden change 
of the field value. 

In order to calculate the change the field needs to real- 
ize the jump, we choose, without losing generality, a > 
so that the region of A > is the basin of the de-Sitter- 
like dominant attractor (c). To meet the constraints from 
nucleosynthesis and structure formation, we require that 
the Quintessential field have well scaled with the back- 
ground before nucleosynthesis. So we choose > 20 
for {VL^{T ~ IMeV) < 0.2). At some recent time just be- 
fore the jump the field (f> = <j)j = ipj and then (j) rapidly 
changes from (pj to (pj + (50 (or from (fj to fjifj)- For 
jumping from the basin of the scaling attractor (6) to 
that of the de-Sitter-like dominant attractor (c) happen- 
ing, we have 



> 



— ln(— 77) — 0, = — In - — — - 

a ^ ' ^ a Xj + af3 



where A, 



fjV] > ln(-?7) ■ 

a 



(30) 



(31) 



Now we shall justify the introduction of the field whose 
value has a sudden change. Below, we suggest a possi- 
ble mechanism, which resorts to spontaneous symmetry 
breaking, to realize the sudden change of the field value. 
We first note that scalar fields are ubiquitous in super- 
symmetric theory of particle physics. Thus it is reason- 
able to assume that a few of them are relevant to the 
cosmic evolution. Considering (p as an effective field, we 
involve three fields with the Lagrangian 

= --f{cT)g^''d^,ipd^ip-- ^° ,, 

2 ^ ^ (ry + e-"v/('T))/3 

-\9''''d^ad^a - ^9^" 8^68,6 - V{a,e) , (32) 
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where 

9) = K - + Ja^^ + + \x'a^e' 

= ^X{M^-a^f + ^m^e^ + ^X'a^9\ (33) 
We conceive /(ct) as 

so that if decouples from a and 6 except at the points 
a = dzcTg (they surely couple to each other through Fried- 
mann equation; nevertheless, when the radiation or mat- 
ter dominates the universe, the coupling through gravity 
is neglectable.). Since /(cr) is more reasonable to be a 
continuous function, a better choice of f{a) might be 

fi'^) = - tanh[a(a2 _ ^2)] ^ ^ » ^ (35) 

We note that V{a;d) is famous for its realization of 
Hybrid inflation models [3l|. In these models, first, a is 
held at the origin, with slow rolling down the potential, 
giving the inflation, and then, when 6 rolls down a critical 
value 9c, (J is destabilized and quickly rolls down from 
to ±M, ending the inflation. Comparing two ways of 
writing the potential, we obtain 

ml = AM^, (36) 

K)' - Jam*. (37) 

And the critical value of 9 is 

9l^ml/X' ^XMyX' . (38) 

For this potential to be viable for current purpose, 9 does 
not necessarily slow roll. Yet we do need a quickly roll 
down from to M after 9 rolls down 9c, which implies 

ml » , (39) 



and we require < ct^ < so that when cr rolls down 
from to ±M, /(cr) changes from 1 to fj. 

Note that when the field (j) jumps to the de-Sitter-like 
regime, V{(f> — iff (a)) will increase and the kinetic term 
of (f will also change. At the same time, V{a, 9) should 
decrease so as to vanish when Quintessence begins to 
dominate the universe. For this scenario to be viable, 
we require that the decrease of V{a, 9) be larger than the 
increase of V^((/> — (pf{a)) and the kinetic term (note that 
when \fj\ < 1, the kinetic term will decrease; however, 
it is easy to show that the decline of the kinetic term 
in this case is small, compared to the increase of the 
V{<p — ff{^)) )- One may worry that this might spoil 
the analysis of the dynamics of (j) above, as this require 
the energy associated with cr and 9 to be comparable with 
that associated with Lp around the jump point. However, 
we argue that it will not, because Lp almost decouples 
to cr and 0, and the energy associated with a and 9 is 
only comparable with that associated with tp when the 
radiation or matter dominates the universe and vanishes 
when the dark energy begins to dominate the universe. 

In summary, we suggest a new approach to construct 
Quintessential dark energy models, with which we first 
propose a relation F = /(A) between F = VV" / {V'Y and 
A ~ —Y'lV, and then figure out the potential V{<^^- It is 
showed that a tracker solution that is different from those 
discovered before can be easily obtained and a solution of 
multiple attractors is also found straightforwardly. Then 
we suggest a scenario that the initial conditions of the 
cosmic scalar field are in the basin of a scaling attractor 
and first, the field evolves toward the scaling attractor 
and then recently, the field jumps out to the basin of a 
de-Sitter-like dominant attractor, giving rise to current 
acceleration. For this scenario to be realized, we intro- 
duce a field whose value changes a certain amount in a 
short time. Then we calculate the change the field needs 
for a successful jump and invoke a mechanism that is 
similar to the case of Hybrid inflation to justify the in- 
troduction of this kind of field. 

We thank Yun-Song Piao, Yi Wang, Zuo-Tang Liang, 
Jian-Hua Gao and Ye Chen for useful discussions. 
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